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Abstract 

The method for the evaluation of the non-resonant thermonuclear function in the Maxwell- 
Boltzmann case with depleted tail is discussed. Closed forms of the analytical results are 
obtained in computational format, and written in terms of the .ff-function in two variables. 
The standard non-resonant cases are extended to Tsallis reaction rates through the pathway 
model. Behavior of the depleted non-resonant thermonuclear function is studied. A compar- 
ison of the Maxwell-Boltzmann energy distribution with a more general energy distribution 
called pathway energy distribution is also done. 

Keywords: Thermonuclear function, pathway model, reaction rate probability integral, H- 
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1 Introduction 

Thermonuclear reactions play a vital role in the theory of fusion physics. The evaluation of 
the reaction rates for low-energy non-resonant thermonuclear reactions in the non-degenerate 
case is done using the principles of nuclear physics and kinetic theory of gases [13] . A nuclear 
reaction in which a particle of type i strikes a particle of type j producing a nucleus p and 
a new particle q is symbolically represented as i + j — > p + q. If and rij are the number 
densities of particles i and j respectively and if the reaction cross section is denoted by cr(y) 
where v is the relative velocity of the particle and f(v) is the normalized velocity distribution, 
then the thermonuclear reaction rate is obtained by averaging the reaction cross section 
over the normalized distribution function of the relative velocity of the particles given by 
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rij=mnj va(v)f(v)dv = n i n j {av) ij . (1.1) 
Jo 

The bracketed quantity {crv)ij is the probability per unit time that two particles of type % and 
j confined to a unit volume will react with each other. For a non-relativistic, non-degenerate 
plasma of nuclei in thermodynamic equilibrium, the particles in the plasma possess a classical 
Maxwell-Boltzmann velocity distribution given by |16j 

SmbdWv = (^) 2 exp 4™ 2 d V , (1.2) 

where /i is the reduced mass of the particles given by fi = ^p^-, T is the temperature, k is 
the Boltzmann constant. Writing in terms of the relative kinetic energy E = we get the 
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Maxwell-Boltzmann energy distribution as (3j [15] 
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Using (1.1) and (1.3) we have, 
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For a non-resonant nuclear reactions between two nuclei of charges Zi and Zj colliding at low 
energies below the Coulomb barrier, the reaction cross section has the form (3j E] 
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where e is the quantum of electric charge, /*& is the Plank's quantum of action and S(E) is the 
cross section factor which is often found to be constant or a slowly varying function of energy 
over a limited range of energy given by |21| [9] 
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Substituting (1.5) and (1.6) in (1.4) we get 
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Putting y = -g? and x = 27r (3^) — ^ 
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Thus the reaction rate probability integral in the Maxwell-Boltzmann case is given by 
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Let us consider a general form of the integral as 



y i/-i e -*v-*v p dy, 7 G C, z > 0, x > 0, p G 



(1.10) 



Physical situations different from the ideal non-resonant Maxwell-Boltzmann case can be 
obtained by modification of the cross section o~(E) for the reacting particles and/ or by the 
modification of their energy distribution. Some of the non standard physical situations are as 
follows 01251 El : 

Non-resonant case with high energy cut-off 
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If the thermonuclear fusion plasma is not in a thermodynamic equilibrium then there is 
a cut-off in the high energy tail of the Maxwell-Boltmann distribution function, then the 
thermonuclear function to be evaluated takes the form 

Il(u,l,x, h= [* y»e-y~ x y~ h dy, x > 0, d < oo. (1.11) 
1 Jo 

The general form of the integral in this case can be taken as 

rd 

-l,z,x,p) = y-r-i e -*v-xv p dy, 7 G C, z > 0, x > 0, d < oo. (1.12) 
Jo 

Non-resonant case with depleted tail 

If we consider an ad hoc modification of the Mawell-Boltzmann distribution which looks 
like a depletion of the tail of the Maxwell-Boltzmann distribution as suggested by Eder and 
Motz [4|, Clayton et al. [2] and Mathai and Haubold [21j which is given by 

1 f°° f, - 1 

I 3 (v,l,l,6,x,-)= / y" e -y-*-*» 2 dy, x > 0, 5 E M+. (1.13) 
z Jo 

We will consider here the general integral of the type 

/ 3 (7 " 1, t, z, 5, x,p)= y-r-ie-ty-zy ^'"dy, (1.14) 

Jo 

where 7 e C, t > 0, z > 0, x > 0, p G M + , <5 G M + . 
Non-resonant case with screening 

The electron screening effects for the reacting particles can modify the cross section of the 
reaction. The reaction rate probability integral in this case will take the form 

1 f°° - - 

h(u,l,b,t,-)= aTe-v-Av+t) 2 dy, x>0,t>0 (1.15) 
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where t is the electron screening parameter. Here we consider the general integral as 

;>oo 

h(j-l,z,x,t,p)= y*- 1 e -* v - x (v+ t )- p dy, 7 G C, z > 0, x > 0, t > 0,p G M+. (1.16) 

The evaluation of the integrals ii, 1^,13 and ^4 in the physical and Astrophysical literature 
are by approximating the integrals by means of the method of steepest descent [3 El [8] . The 
closed forms of the integrals I±, 1^,13 and I4 in terms of Fox's H- function and Meijer's G- 
function can be seen in a series of papers by Mathai and Haubold, see for example Haubold 
and Mathai[16J, Mathai and Haubold |17l I22| etc. In the present paper we will consider the 
integral -Z3 in the depleted case in detail and obtain the closed form evaluation of the function 
by a different method. Also we extend the integral to a more general case than the Maxwell- 
Boltzmann case using the pathway model introduced by Mathai in 2005. 

The paper is organized as follows: In the next section we consider the general form 
of the non-resonant reaction rate probability integral in the Maxwell-Boltzmann case with 
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depleted tail and obtain the closed form via the ii-function in two variables. A more general 
form of the depleted non-resonant thermonuclear function is obtained by using the pathway 
model in section 3. Section 4 is devoted to study the behaviour of the depleted non-resonant 
thermonuclear function in the Maxwell-Boltzmann and Tsallis case and compare the Maxwell- 
Boltzmann energy distribution with a more general energy distribution. Concluding remarks 
are included in section 5. 

2 Closed form representation of the standard non-resonant 
thermonuclear function with depleted tail 

In this section we evaluate the integral ^3(7 — 1, t, z, S, x, p) and give a representation for 
it in terms of H- function in two variables. 

For non-negative integers mi,m2,m3,ni,n2,n3,pi,p2,P3,qi,q2,Q3 such that < mi < 
c/1,0 < 77J2 < <Z2>0 < 7713 < 53,0 < Ti2 < P2,0 < 773 < ps, for ai,bj,Cj,dj,ej, fj £ C and for 
ctj, Pj, Aj, Bj,Cj, Dj, Ej, Fj e M + = (0, 00), the //-function in two variables is defined via a 
double Mellin-Barnes type integral in the form 
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and x and y are not equal to zero, and an empty product is interpreted as unity. The contour 
L\ is in the si-plane which runs from b\ — ioo to <5i + too, which separates all the poles of 
T(bj+ f5jS\ + BjS2) and T(dj+DjS\) to the left and all the poles of T(l — Cj — CjS\) to the right. 
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The contour L 2 is in the S2-plane which runs from <5 2 — ioo to 62 + ioo, which separates all the 
poles of T(bj + fijS\ + B3S2) and + FjS2) to the left and all the poles of T(l — ej — EjS2) 
to the right. The H- function in two variable given in (2.1) will have meaning even if some 
of these quantities are zeros. For details about the contours and existence conditions see 
Srivastava et al.[28j, Mathai and Saxena |27j . The details of the //-function and G- function 
in one variable can be seen in |19| 1251 [26] 

Let the function /(xi,x 2 ) be defined in = (0, +00) x (0, +00). Then the Mellin 
transform of a function f(xi,x 2 ) in points (si,s 2 ) £ C 2 is defined as 
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with the inverse 



f(xi,x 2 ) = . .- 2 / / Mf(si,S2)x^ Sl X2 S2 dsids2- (2.6) 

{Z7TI) JSj-ioo J 62-100 

The conditions under which the (2.5) and (2.6) are valid have been discussed by Fox [10J and 
Hai and Yakubovich |12|. 



Now consider the integral ^3(7 — 1, t, z, 6, x, p) given in (1.14). We evaluate this integral 



by using the Mellin transform technique for two variables. Using (|2.5|) and 
f{t,z) = 13(7- l,t,z,6,x,p) = 

we have, 
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Changing the order of integration due to the uniform convergence of the integral, we get 
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Putting xy p = u we get, 
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Taking the inverse Mellin transform using (2.6) we obtain, 
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where flo'i-'o'l-'o'l * s an -f^- function in two variables defined as in (2.1 ). If ^ is an integer then 
put -= m,m = 1,2, ■■■ . Then using the multiplication formula for gamma function defined 
by [MIES] 

T(mz) = (2vr) i ^ 1 m m2 -5r(z)r (z+-Y--v(z+ ^— , (2.10) 

V ml \ ml 



where z£C,z/0, —1, —2, ... and m a positive integer, we have (3.17) as 
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(2.11) 



For the non- resonant case with depleted tail we have j = 1 + u, p = \, Then by using the 
duplication formula for gamma functions we obtain 
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Next we obtain the extension of these results by using the pathway model of Mathai. 



3 Extension of the non-resonant thermonuclear function with 
depleted tail through pathway model 

In this section we try to extend the non-resonant reaction rate probability integrals to a 
more general case. The extension is done by using the pathway model introduced by Mathai 
in 2005 [20^23] . This model was first introduced for the matrix variate case but here we make 
use of the scalar case of the model for extension of the results. By the pathway model one 
can move between three different functional forms namely the generalized type-1 beta form, 
generalized type-2 beta form and the generalized gamma form. The pathway model for the 
real scalar case is defined as follows: The generalized type-1 beta form of the pathway model 
is given by 

h(x) = c i x' y ~ 1 [l -a{l-a)x s \ x=s, a > 0, 5 > 0, 1 - o(l - a)x s > 0, 7 > 0, a < 1 (3.1) 

where a is the pathway parameter. This is the case of right tail cut-off. For a = 1,7 = 1,6 = 1 
we get the Tsallis Statistics for a < 1 pH EH [30] . For a > 1 

f 2 (x) = c 2 x 7 - 1 [l + a(a - l)x s ]~^, < x < 00 (3.2) 

is a generalized type-2 beta form of the pathway model. Here also for 7 = l,a = l,<5 = l 
we get the Tsallis Statistics for a > 1 [HI [291 SOI- Superstatistics of Beck and Cohen [1] is 
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obtained for a = 1, 5 = 1. As a — > 1 the functions given in (3.1) and (3.2) will reduce to the 
generalized gamma form of the model given by 



fs(x) = cgsT-V* 8 ,*>0. 



(3.3) 



Here ci, C2 and C3 are the normalizing constants if we consider the above functions as statistical 
densities. Many statistical densities come as particular cases of the above three functional 
forms, see Mathai[20| and Mathai and Haubold |23} I24j for details. By using the principles of 
pathway model we can obtain a new energy distribution given by 
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for a > 1, — ^-j- — I > 0, which is more general than the Maxwell-Boltzmann energy distri- 
bution defined in (1.3). As a — > 1 we obtain the Maxwell-Boltzmann energy distribution. 
Substituting the pathway distribution (3.4) in (1.1) and using (1.5) and (1.6) we obtain the 
reaction rate probability integral in the extended form denoted by fjj as 
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This is the extended non-resonant thermonuclear function in the Maxwell-Boltzmannian form. 

1 2 

Putting y = -§f and x = 2tt (^w) 5 ^ e , we obtain the above integral in a more simplified 
form as 
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for a > 1, — — I > 0. The integral to be evaluated in this case is of the form 

1, 



I\ a {v, 1,X, 



y"[l + {a-l)y]- 



dy. 



(3.7) 



A more general integral to be evaluated in the extended Maxwell-Boltzmann form can be 
taken as 



ha{l -1,2 
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Other general integrals to be evaluated are 

,-d 
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which are the extended cut-off case, extended depleted case and extended screened case re- 
spectively. Among these integrals the closed form representations of hail — l,z,x,p) an d 
I2a.il ~ ^-i z i x iP) m terms of Fox's H- function can be obtained as in [144 115] 
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For the case of astrophysical interest, the extended Maxwell-Boltzmann case or the Tsallis 
reaction rate can be obtained as 
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and the extended cut-off case can be obtained as 
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where G^n is the Meijer's G-function, see Mathai |19| . Mathai and Saxena [25] or Mathai and 
Haubold [25] for details. The detailed evaluation of the integrals in terms of //-function and 
their special cases in Meijer's G-functions can be seen in Haubold and Kumar |14[ll5j. Kumar 
and Haubold [E]. The integral hail— 1, z,x,t, p) can be obtained in terms of hail~ 1> z -,x,p) 
and I^il — 1, z, x, p) by some basic arithmetic procedure. 
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Here we will evaluate the integral /3a (7 — 1, t, z, 5, x, p) and obtain the closed form repre- 
sentation in terms of ^/-function in two variables. For, let us consider the integral 



g(t,z) =I 3a = I yi-\l + (a - l^-^e-^-^dy. 
J 

We will evaluate this integral also by using the Mellin transform technique as in the case 
discussed in the previous section. We have, 
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Jo Jo Jo 

Changing the order of integration and simplifying using suitable substitution we get get 
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where > 0,3i(s 2 ) > 0,3ft — si) > 0. Then simplifying exactly as in the previous 

case we get, 
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where ^o'l l'l-o'i ^ s an ^-function in two variables defined as in (2.1 ). If - = to, to = 1, 2, • • • 
then by using ( |2.10[ ) we get 
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For the extended non-resonant case with depleted tail we have 7=1 + 1/, p= \, we have, 
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In the next section, we compare the standard non-resonant thermonuclear function with 
depleted tail with the extended non-resonant thermonuclear function with depleted tail. 

4 Comparison of the extended results with the standard re- 
sults 

Here we try to compare the results obtained in the standard and extended non-resonant 
thermonuclear functions in the standard and extended case. In the Mellin-Barnes integral 
representation of (3.19) given by 



h a {v, 1, l,6,x, 
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if we take the limit as a 

mm 



1 then by using the asymptotic expansion of gamma function 



F(z + a) ~ (2vr)Sz 2 + a -i e 



oo, \arg(z + a)\ < ir — e, e > 0, 



(4.2) 



where the symbol ~ means asymptotically equivalent to, we get (2.12). Next we compare the 
Maxwell-Boltzmann energy distribution with the Pathway energy distribution. Figure l.(a) 
shows the Maxwell-Boltzmann energy distribution for the value of kT = 100, 200, 300. As 
we increase the value of kT it is observed that the function is heavy tailed and less peaked. 
Figures l.(b),(c) and (d) show the pathway distribution for kT = 100,200,300 respectively. 
Jpd(E) is plotted for a = 1, a = 1.1, a = 1.2, a = 1.3, a = 1.5 and a = 1.6. 
SI 

° % W=100 




jtr=30o 




(b) 
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Figure l.(a) f M BD(E) for kT = 100,200,300. 

(b) fpo{E) for kT = 100, a = 1, a = 1.1, a = 1.2, a = 1.3, a = 1.5 and a = 1.6 

(c) Jpd{E) for kT = 200, a = 1, a = 1.1, a = 1.2, a = 1.3, a = 1.5 and a = 1.6 

(d) fpD(E) for A;T = 300, a = 1, a = 1.1, a = 1.2, a = 1.3, a = 1.5 and a = 1.6 

From the graphs it can be observed that the pathway energy distribution {fpp>{E)) is more 
general than the Maxwell-Boltzmann energy distribution (fMBD(E)). We can retrieve the 
Maxwell-Boltzmann energy distribution from pathway distribution as a — )■ 1. As we increase 
the value of kT in fpp>(E) we observe that the function becomes thinker tailed and the peak 
is reduced. 



5 Conclusion 

An attempt has been made to change the energy distribution of the ions in the plasma 
from the Maxwell-Boltzmann case. By this change of using the pathway energy distribution, 
more unstable and chaotic situations are covered whereas the standard Maxwell-Boltzmann 
situation is retrieved by letting a — > 1. It may be noted that even a small deviation of the 
energy distribution with a produce a dramatic effects on those nuclear reaction rates whose 
main contribution comes from the high energy tail of the distribution which can be observed 
from the Figure. The standard and extended non-resonant thermonuclear functions with 
depleted tail are evaluated by using the Mellin transform technique helped to obtain more 
convenient closed form representations. The figures are plotted by using Maple 14 under 
Microsoft Windows XP platform. 
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